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ABSTRACT 

We used the VLA to search for ^He^ emission from two Galactic planetary 
nebulae (PNe): NGC6572 and J 320. Standard stellar models predict that the 
^He/H abundance ratios for PNe should be 1-2 orders of magnitude higher than 
the primordial value (^He/H ~ 10~^ by number) determined from Galactic H ii 
region abundances and confirmed by WMAP cosmic microwave background re- 
sults. Chemical evolution models suggest that fewer than 5% of all PNe enrich 
the interstellar medium (ISM) with ^He at the level of standard stellar mod- 
els. Our target PNe are therefore anomalous in that they were selected from 
a sample deliberately biased to contain objects with properties that maximized 
the likelihood of a ^He detection by the VLA. We have detected the 8.665 GHz 
hyperfine ^He^ transition in J 320 at the A a level. The ^He/H abundance ratio 
is 1.9 X 10~^ with roughly a factor of two uncertainty. For NGC6572 we find an 
upper limit of ^He/H < 10~^ . This detection of ^He in J 320 makes it the second 
PN known to have an anomalously high ^He abundance confirming that at least 
some low-mass stars produce significant amounts of ^He that survives to the PN 
stage and enriches the ISM. 

Subject headings: ISM: abundances — planetary nebulae: general — radio lines: 
ISM 



1. Introduction 

The light elements, (D), ^He, ^He and ''Li, are expected to be produced during the 
era of primordial nucleosynthesis and then further processed by stars over many generations 
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(Reeves 1974; Boesgaard & Steigman 1985). Observations of the light elements can in 
principle be used to constrain models of stellar and Galactic evolution and cosmology. For 
each light element the quantity of interest is the abundance ratio relative to hydrogen (e.g., 
^He/H). In some cases converting the observables into an abundance ratio can be as difficult 
as making the observations. 

For many years the light elements have been used to constrain standard Big Bang nucle- 
osynthesis (SEEN) models (Copi et al. 1995). SEEN contains only one free parameter, the 
baryon-to-photon ratio, rj\^. Observations of optical recombination lines of ^He in metal poor 
galaxies were first used to probe SEEN (Peimbert & Torrcs-Pcimbcrt 1974). Observations 
of the resonance line of ^Li in metal poor Halo stars revealed a constant ^Li/H abundance 
ratio, called the "Spite Plateau", that was interpreted as the primordial abundance and 
used to determine values of 77b (Spite & Spite 1982). The advent of large aperture, optical 
telescopes provided the sensitivity to detect the Lyman-series lines of D in metal poor QSO 
absorption systems (Tytler et al. 1996). More recently the ^Hc"'" hypcrfinc transition was 
used to estimate 77b by measuring ^He/H in H 11 regions throughout the Galaxy. A "^He 
Plateau" was discovered and from it a primordial abundance inferred (Eania et al. 2002, 
and shown in Figure 1). Since the primordial D/H abundance ratio is more sensitive to ?7b, 
and the processing of D in stars is thought to be understood, it became the baryometer of 
choice by many investigators (Schramm & Turner 1998, for example). Nonetheless, all the 
light element abundances should agree with SBBN within the uncertainties. 

Observations of the cosmic microwave background (GME) provide an independent method 
of determining ?]b- GME anisotropy measurements made with the Wilkinson Microwave 
Anisotropy Probe (WMAP) were used by Spergel et al. (2003) to derive the primordial 
baryon-to-photon ratio. Adopting the CMB results for 77b gives the primordial light element 
abundances (Romano et al. 2003; Cyburt et al. 2003; Coc et al. 2004; Izotov & Thuan 
2004; Steigman 2005). The D and ^He abundances are in good agreement with those derived 
from the WMAP observations. Depending on the adopted uncertainties the measured ^He 
abundance is either in agreement with GME data or slightly lower than expected (Peimbert 
et al. 2003; Olive & Skillman 2004). The ^Li abundance is significantly lower given the low 
dispersion of these measurements. A possible resolution of this result is that the ^Li abun- 
dances are depleted in metal poor Halo stars (Charbonnel & Primas 2005, and references 
within) . 

Figure 1 shows the ^He/H abundance ratio by number as a function of Galactocentric 
radius (Eania et al. 2002). The solid circles (and crosses) are ^He abundances for H 11 
regions located throughout the Galaxy. Since H 11 regions are very young, they measure 
the ^He/H abundance ratio at the present time. Not plotted are the ^He/H abundances for 



-3- 



protosolar material and the local ISM derived from the Galileo Probe (Mahaffy et al. 1998) 
and Ulysses (Gloeckler & Geiss 1998), respectively, that are in accord with the H ii region 
^Hc abundances. The solid triangle is the ^He abundance for the planetary nebula NGC 3242. 
Early stellar and Galactic evolution models predicted that solar type stars would produce 
and eventually expel into the ISM prodigious amounts of ^He, creating an easily measurable 
abundance gradient with Galactocentric radius or metallicity (Rood et al. 1976). The high 
^He abundance in NGC 3242 was consistent with these early models but inconsistent with 
the H II region ^He abundances, stimulating a new generation of models. 

In standard stellar models "^Hc is produced on the main sequence outside the main 
nuclear burning regions where p-p burning terminates in '^Hc. Longer lived low mass stars 
produce more ^He. This ^He is mixed into the convective envelope at the first dredge up on 
the lower red giant branch. Once in the envelope the ^He survives, so any stellar winds and 
the final PN should be substantially enriched in ^He. 

Standard stellar models include only mixing produced by thermal convection. Obser- 
vations have long indicated that extra mixing processes must be taking place, e.g., over- 
shooting beyond convective radiative boundaries or rotationally induced turbulence (Zahn 
1992). Charbonnel et al. (1998) argued that this rotationally induced extra- mixing occurs 
in low mass stars above the luminosity function bump produced when the hydrogen burning 
shell crosses the chemical discontinuity left by the retreating convective envelope early in 
the red giant branch phase. This hypothesis was supported by observed increases in the 
^^C abundance above the luminosity function bump. The same mixing that would increase 
the ^^C abundance would bring the '^He down to values comparable to the primordial value. 
Further extra-mixing could take place on the second ascent asymptotic giant branch. In 
stars which have undergone extra mixing the material expelled into the ISM during the PN 
stage is therefore not enriched in ^He as originally thought. Charbonnel & do Nascimento 
(1998) predict that 96% of all low-mass stars will undergo this mixing process. The expected 
^He yields coupled with Galactic evolution models (Tosi 2000; Chiappini et al. 2002; Galli 
2005) are consistent with the "^He Plateau" shown in Figure 1 and described by Bania et 
al. (2002). 

The horizontal grey band in Figure 1 is the ^He/H primordial abundance ratio predicted 
by SBBN calculations (Buries et al. 2001) using the WMAP 7/b value. The CMB abundance 
lies at the lower envelope of the ^He/H abundances. This and the absence of an observable 
negative gradient in the ^He abundances with galactocentric radius led Bania et al. (2002) 
to argue that the stellar contribution to ^He is positive definite but very small. The scatter 
in the H ii region abundances could be real (e.g., azimuthal abundance fluctuations). But 
this scatter could also be measurement error (e.g., the ionization correction could be un- 
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derestimated for some H ii regions) . The scatter presumedly obscures a yet to be observed 
gradient. Given this, a hmit on BBNS production can be determined from the best deter- 
mined abundance in the outer Galaxy. The rightmost solid point in Figure 1, the H ii region 
S209, fits this requirement, and Bania et al. (2002) took the S209 abundance as their best 
estimate of the primordial abundance. This value is almost identical to the WMAP value. 

There arc two observational paths to improving our understanding of the evolution of 
^He. The first is to reduce the scatter in the H ii region abundances by observing a sample 
chosen because of their simple physical structure and optimum location in the Galaxy. The 
second is to observe the stellar sources of "^He and directly measure the "^He/H abundance 
in PNe. 

Previous observations of six PNe were made using the Max Planck Institut fiir Radioas- 
tronomie 100 m telescope (Balser et al. 1997). The PN sample was chosen based on criteria 
that would maximize detection, including diagnostics that indicated less mixing. ^He^ was 
only detected in NGC 3242, shown as the sohd triangle in Figure 1. It therefore appears that 
in at least one object significant amounts of ^He in a low-mass star survives until the PN 
phase, and that PN and its earlier stellar winds have enriched ^He in the ISM. We now seek 
to understand why some stars produce ^He and others do not. 

Observations of ^He"*" in PNe are challenging because of the weak intensity of the ^He"*" 
emission line coupled with their small angular size. We are increasing the size of our PN 
sample using the Very Large Array (VLA), the 100 m Green Bank Telescope (GBT), and 
the 305 m Arecibo telescope. Our target sample is purposely biased to maximize detections. 
Our target selection criteria have evolved with time. Currently they are: 

• Any mixing in the progenitor star beyond that included in standard stellar models is 
likely to destroy ^He. Indications of mixing include: high ^^C/^^C, N/0, and ^He/H 
abundances. To be included as one of our targets a PN must have a ^He/H abundance 
ratio < 0.125 by number and [N/0] < —0.3 with a value < —0.6 considered even 
better. Unfortunately the ^^C/^^C abundance ratio is known for only a small sample 
of PNe (Clegg et al. 1997; Palla et aL 2000, 2002; Balser et al. 2002). 

• The ^He abundance in the progenitor grows as main sequence lifetime increases, hence 
PNe from older stellar populations are more promising candidates. The Peimbert Class 
(Peimbert 1978; Faiindez-Abans & Maciel 1987) gives an estimate of the progenitor 
age (as well as abundance criteria similar to ours). Our targets are chosen from classes 
lib. III, and IV. 

• Since we observe ^He'*', our target PNe must have excitations such that most of the 
He is singly ionized. 
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Here we discuss ^He+ VLA observations of two PNe: NGC 6572 and J 320. Both objects 
are expected to have high ^He/H abundance ratios; they were selected using criteria similar 
to those above. In addition, the angular sizes of these PNe are well matched to the VLA 
synthesized beam at 8.7 GHz. 

2. Observations 

Table 1 summarizes the VLA^ D-array configuration observations. NGC 6572 was ob- 
served on five different days. The correlator spectrometer was setup to observe the ^He"*" 
hyperfine transition at two orthogonal polarizations at a rest frequency of 8665.65 MHz. The 
data were Hanning smoothed on-line. During part of the first day the H91a transition was 
observed at a rest frequency of 8584.81MHz. J 320 was observed on seven different days. 
The correlator was setup to observe the •^He^ and HQlo; transitions simultaneously at two 
orthogonal polarizations each. The data were Hanning smoothed off-line. 

The data were edited and calibrated using the Astronomical Image Processing System 
(AIPS). The calibrator 3C286 was used to calibrate the intensity scale. Because the '^He''' 
emission line intensity was expected to be weak the bandpass had to be carefully calibrated. 
The spectral baselines are typically much flatter for interferometers than single-dish tele- 
scopes since any telescope dependent effect will be removed when the data from different 
antennas are cross correlated. But at the sensitivity levels required for our observations the 
VLA bandpass is not stable with time; amplitude variations as large as 1% have been mea- 
sured with an approximate width of 3 MHz.^ These baseline ripples are caused by reflections 
in the waveguides and move across the bandpass at about 0.4 MHz per hour. We therefore 
observed a bandpass calibration source approximately every 15 minutes for a duration of 
about 5 minutes. A bandpass model was calculated by interpolating the bandpass solutions 
over time. We are confldent that we have removed most of this structure since the measured 
spectral rms noise is consistent with the expected thermal noise. 

The continuum intensities of NGC 6572 and J 320 were sufficient to perform the standard 
self-calibration technique. We applied the self-calibration algorithm to each day's observa- 
tions after the standard phase and bandpass calibration. This procedure typically improved 
the signal-to-noise ratio by a factor of five in the continuum. The task DBCON was used to 
average the data. The self-cahbration technique was applied to the flnal average. The task 
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UVLSF was used to subtract the continuum emission in the {u, T;)-plane using the hne-free 
channels. This fits a hnear basehne to the specified channels in the real and imaginary parts 
of the visibilities. Continuum and line images were produced using the CLEAN algorithm. 
Table 1 lists the rms noise in these images. The continuum images arc completely limited 
by the finite dynamic range of the strong continuum and thus are not limited by thermal 
fiuctuations. The line data, after subtraction of the continuum, are limited by thermal noise. 
Thus the continuum rms values are larger than the line data. 



3. Observational Results 

The continuum images of NGC 6572 and J 320 are shown in Figures 2 and 3, respectively. 
The continuum emission is well modeled by a Gaussian function. The task JMFIT was used 
to fit a single- component, two-dimensional Gaussian model to the continuum image for each 
PN. Table 2 summarizes the results. Listed are the source name, the B1950 equatorial 
coordinates of the peak radio emission, the peak flux density, the integrated flux density, 
and the angular size. 

We analyzed the H91q; and "^Hc^ spectral line data cubes using four different methods: 
(1) the ISPEC task was used to produce spectra toward the continuum peak emission within 
a single resolution element; (2) averages over different spatial areas were made in order to 
increase the signal-to-noise ratio; (3) the images were convolved with Gaussian functions that 
had half-power beam-widths (HPBWs) larger than the synthesized beam; and (4) the task 
IRING was used to compute averages over concentric annuli to produce a cumulative flux 
density versus radius relation. This latter method averages over both position and velocity. 

We detected HQlo; line emission in both PN; only J 320 has ^He"*" emission. Spectra 
are shown in Figures 4 and 5 for NGC 6572 and in Figures 6 and 7 for J 320. The solid line 
histograms are the data and the dashed lines are models described in §4. The additional 
spectra shown for J 320 were made by convolving the H91q; data cube with a HPBW of 20" 
(Figure 6) and the ^He^ data cube with HPBW's of 20, 30, and 40" (Figure 7). Comparing 
the ^He'*' line intensity to the noise level in the line free region in Figure 7 shows that the 
detection has a signal to noise ratio of S/N ~ 4. Figure 8 plots the cumulative ^He"*" flux 
density from the ^He'*' channels (12-22) versus radius for J 320. The ^He"*" flux density 
gradually increases until a radius of 15-20". For larger radii the flux drops because of 
missing short spacing data. In the configuration used for these observations the VLA has 
no sensitivity to ^Hc'^ emission with radii larger than about 30". The plot also shows the 
cumulative fiux density for the channels not containing the line emission (6-11 & 23-26). A 
comparison of the line and line- free results gives a 5'/A'" of ~ 9 for the ^He^ detection. This 
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value is larger than that obtained from the spectra because the ^He"*" emission is averaged 
over both space and frequency. The detected continuum and recombination hue emission 
are both confined to a smaller area within a radius of 5-10". Apparently much of the ^Hc^ 
emission from J 320 emanates from a large, low density halo. The model shown in Figure 8 
has a 35" halo (see §4.2). This halo model fits the data well except for larger radii where it 
is not affected by the missing short baselines. 



4. The ^He Abundance 

Estimates for the ^He/H abundance ratios in the PNe were made by numerical modeling. 
We used the radiative transfer code NEBULA to calculate synthetic spectra of the radio con- 
tinuum, radio recombination line, and ^He"^ line emission from a model nebula (Balser et al. 
1999a). The model is comprised solely of ionized hydrogen and helium gas. A density and 
ionization structure is specified within a three-dimensional Cartesian coordinate system. For 
each numerical cell the following physical properties are specified: electron temperature (Tg), 
electron density (ue), and hehum ionization (^He''"/H+, ^He''""'"/H+, ^He''"/H+). The mod- 
eled nebula consists of homogeneous concentric shells which have a radial expansion velocity 
(Balser et al. 1997, 1999b). The line profile is therefore broadened by three components: 
thermal, microturbulent, and large-scale expansion. The thermal motions are determined 
from the electron temperature, whereas the microturbulent motions, by definition smaller 
than the numerical cell, are set by a single free parameter that is constrained by the observed 
linewidths. The ^He"*" and recombination lines are assumed to be in LTE. Non-LTE effects 
and pressure broadening from electron impacts are negligible for these PN models. Syn- 
thetic spectra are generated by calculating the radiative transfer for the line and continuum 
emission as a function of frequency from the back to the front of the numerical grid. The 
NEBULA outputs are the two-dimensional continuum brightness distribution and spectral 
data cubes for ^He'*' and the recombination lines. These model data are then convolved with 
a specified telescope beam so they can be compared directly with observations. 

We use an iterative procedure to derive an estimate for the ^He/H abundance. First 
the physical properties of a source arc constrained by the continuum and recombination line 
emission as well as data from the literature. A fiducial ^Hc/H abundance is then chosen and 
synthetic spectra are calculated with NEBULA. The "^He/H abundance is iterated until the 
model spectra match the observations. A grid of models is used to assess the uncertainty in 
the derived '^Hc/H abundance. 

The planetary nebulae properties used as input for NEBULA are summarized in Table 3. 
The PNe distance, expansion velocity, and helium ionization structure are taken from the 
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literature (see §4.1-4.2). The VLA data are used to model the PNe as single, homogeneous, 
spheres (Balser et al. 1995, 1999a). From these models we determine the spherical angular 
size, electron temperature, and electron density. Below we discuss each PN separately. 

4.1. NGC6572 

NGC6572, discovered by Struve in 1825, is a well studied planetary nebula (Acker 
et al. 1992). Unhke most PNe the distance to NGC6572 has been accurately measured 
by combining the angular expansion rate, determined from high spatial resolution radio 
continuum images at two or more epochs, with the Doppler expansion velocity (Masson 
1986). We adopt a distance to NGC6572 of 1.2 kpc (Kawamura & Masson 1996). The 
expansion velocity is directly measured from spectral line data that resolves the PN shell 
(Masson 1989; Acker et al. 1992). The helium ionization properties are taken from Cahn 
et al. (1992). 

The VLA data were used to model a single, homogeneous sphere with an angular size 
of 7'.'9. Although NGC6572 is just resolved by our VLA data this size is consistent with 
other estimates (Acker et al. 1992). NGC6572 has been imaged at higher resolution in 
the radio (Masson 1989), near infrared (Hora et al. 1990), and optical (Miranda et al. 
1999). Morphologically, NGC6572 has a double-lobed, elliptical structure with a central 
minimum. We assume an inner shell size of 2". The electron temperature and density are 
found to be 10,300 K and 22,000 cm~^, respectively. Various values for Tg and rig can be 
found in the literature. These values depend on the spectral lines employed, the method 
used to determine these physical parameters, and the volume within the nebula sampled 
since real temperature and density fluctuations are present (Mathis et al. 1998). For 
example, electron temperatures determined from optical forbidden line ratios are weighted 
toward higher temperature regions whereas radio recombination line-to-continuum ratios are 
weighted toward lower temperature regions (Mansfield 1969). Nevertheless, our estimates 
are in good agreement with Zhang et al. (2004), who use the optical hydrogen recombination 
spectrum to determine and Ue- 

The model results are shown as dashed curves in Figures 4 and 5 for the H91q; and 
^He"*" lines, respectively. The H 91q; model profile is a good fit to the observations. The data 
are consistent with an upper limit for the ^He/H abundance ratio of < 10~^. 
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4.2. J 320 

Discovered in 1913 (Jonckheere 1913), J 320 has a radio continuum flux density about 
50 times weaker than NGC6572; thus it is not as well studied. The morphology and 
kinematics of J 320 have recently been investigated using long slit spectroscopy with the 
Anglo- Australian Telescope (AAT) and H a. images with the Hubble Space Telescope (HST) 
(Harman et al. 2004). The central nebula consists of two (possibly three) bipolar lobes 
expanding at ~ 46 km sec~^ with different orientations. The morphology is thus better clas- 
sified as a poly-polar PN. Surrounding the central nebula are high speed knots that may not 
have played a significant role in forming the central nebula. 

The distance estimates to J 320 vary from 2 to 6 kpc using different methods (Harman et 
al. 2004, and references within). Here we adopt a distance of 5 kpc. The hehum ionization 
properties are taken from Cahn et al. (1992). The VLA radio continuum and recombination 
line data are used to determine the angular size, electron temperature, and electron density. 
The angular size is 7'.'4, consistent with previous estimates (Acker et al. 1992; Tylenda et 
al. 2003). We assume an inner shell size of 2". The electron temperature and density 
are 7,500 K and 1,500 cm"^, respectively; these values are lower than predicted by optical 
forbidden fine ratios (Barker 1978; Cahn et al. 1992), as might be expected for a nebula 
with temperature and density fluctuations. 

Because the ^He^ emission extends beyond the continuum emission, we infer that a low 
density halo exists for this PN (see Figure 8). Halo components are not uncommon for PNe 
but are difficult to detect in optical (Ha emission) or radio continuum (free-free emission) 
images since these tracers of ionized gas are proportional to the emission measure (oc n^) 
and halo densities are typically low. The ^He"*" emission, however, is sensitive to the total 
column density (oc ^e) . Therefore it is plausible to detect the halo in ^He^ emission but not 
radio continuum emission. The same situation obtained for NGC 3242 (Balser et al. 1997, 
1999b). 

The ^He"*" and radio continuum data together with NEBULA models were used to 
constrain some of the halo properties. To estimate the size of the halo we modeled the 
^He^ emission using NEBULA for different halo sizes. Figure 9 plots the cumulative ^He"*" 
flux density versus radius for the VLA data (solid curve) and for flve different halo models 
(dashed curves). The electron density was assumed to be 100 cm~^ for all halo models and 
the cumulative ^He^ flux density was scaled to match the VLA data since we are only 
interested in the relative shape. Quahtatively a halo size of 35" best fits the data. We set a 
limit for the halo electron density by calculating the continuum emission for different values 
of rif,- Figure 10 plots the cumulative continuum fiux density versus radius for the VLA data 
(solid curve) and for five models (dashed curves) with different halo electron densities and 
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a halo size of 35". For halo densities between < < 50cm~^ the predicted continuum 
emission is within 5% of the measured VLA continuum. The upper limit of Ug = 50 cm~^ 
sets a lower limit to the predicted "^He/H abundance ratio. The expansion velocity of the 
halo is constrained by the ^He^ profile shape. Given the low signal-to-noise ratio of the 
^He^ line and the complex morphology of J 320, it is not clear how to interpret the double 
profile in Figure 7. In our models the Doppler thermal and microturbulent broadening 
together produce a Gaussian line shape. An expansion broadening, however, produces either 
a square-wave profile for an unresolved source or two separate peaks for an optically thin, 
resolved source. Nevertheless, in the context of our spherically symmetric expanding shell 
model we use the overall width of the ^He"*" profile to estimate a halo expansion velocity of 
25 km sec~^. 

Based on these constraints, we adopt a halo size of 35" , an electron density of 50 cm~^ , 
and a uniform expansion velocity of 25kmsec~^. Because we have no information about 
these physical properties in the halo, the electron temperature and ionization properties in 
the halo are assumed to be the same as the core values. Since the ^He/H abundance ratio is 
only weakly dependent on the electron temperature, having to use an assumed value is not 
a serious compromise. Furthermore, since the source of ionization is a hot central star all of 
the helium should be singly ionized. 

The results are summarized in Figures 6-8. The solid lines arc the data whereas the 
dashed curves arc the model results assuming a ^He/H abundance of 1.9 x 10""^. Overall 
the model fits the data reasonably well. The major uncertainties in determining the ^He/H 
abundance ratio in J 320 are: (1) the measured line intensity; (2) the distance; (3) the 
halo electron density; and (4) the halo size. The ^He/H abundance ratio is proportional to 
Tj^RInJ^, where 31 is the ^He'*' fine intensity and ii!sun is the distance (Rood et al. 1984). 
Depending on how the signal-to-noise ratio is estimated, we detect ^He"'" emission at a S/N of 
~ 4 or 9. Taking the more conservative value of 4, the uncertainty in the ^He/H abundance 
is ± 0.5 X 10^'^. An error of 1.5 kpc in the distance gives an error in the abundance of 
± 0.3 X 10^'^. The VLA continuum image limits the halo electron density to less than 
50 cm^"^. Smaller values of He are possible and will increase "^He/H. For example, a halo 
electron density of 10cm~^ will increase ^He/H by 50%. The halo size is constrained by 
the distribution of ^He"*" emission. But a larger halo may exist and not be detected by 
our observations. For example, if the halo size extends beyond about 60" it would not be 
detected in ^He"*" emission with the VLA. A larger halo will reduce our predicted value of 
^He/H. For example, a halo size of 100" will decrease the ^He/H abundance by 50%. None of 
these errors is likely to be normally distributed, so we cannot combine them in a meaningful 
way. Grudely, we estimate that the abundance error is roughly a factor of two. 
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5. Discussion 

The most important result of this paper is that we have detected '^He^ in another 
planetary nebula, J 320, with a "^He/H abundance ratio roughly two orders of magnitude 
larger than the primordial value ('^He/H~10~^ by number) determined from Galactic H ii 
region abundances and confirmed by WMAP cosmic microwave background results. Merely 
detecting ^He"*" in J 320 tells us that its stellar progenitor evolved in a manner well described 
by standard stellar models: it produced ^He during its main sequence lifetime; it did not 
undergo any extra-mixing; it subsequently expelled ^He enriched gas into the ISM. To be 
consistent with chemical evolution models, the existence of the ^He Plateau demands that 
the majority of PNe are not net producers of ^He. J 320 now joins NGC 3242 as only the 
second known example of those rare planetary nebulae which enrich the Galaxy in ^He. 

Our long term goal is to identify some surrogate characteristics which allow us to deter- 
mine what fraction of planetary nebulae enrich the ISM in '^He without having to measure 
^He itself. Now that the sample of '^He PNe has a population of two we can begin this 
attempt. There are some intriguing similarities between NGC 3242 and J 320. They both 
satisfy our selection criteria: each has low He/H and [N/0] abundance which indicates little 
or no extra-mixing; each has most of its He in a singly ionized state; each originates from an 
old stellar population as indicated by its Peimbert Class (NGC 3242 is lib; J 320 is IV). 

Superficially, all Peimbert Class lib. III & IV PNe should be good '^He targets. Class lib 
alone encompasses about 35% of all the known PNe (Quireza, private communication). In 
order to not over-produce ^He, there must be something in addition to either the Peimbert 
Class or our other selection criteria (or both) that determines which PNe have high ^He 
abundances. With the detection of ^He in J 320 we have a first hint as to what that parameter 
might be. Both detected PNe have low density halos which apparently produce the bulk of 
the '^He''' signal that we detect (Balser et al. 1999b, and herein). Both have lineshapes that 
are consistent with optically thin expanding shells. Are low density halos required for high 
^He, or do they simply make it easier to detect if present? 

6. Conclusions 

We detected ^He^ emission with the VLA in the planetary nebulae J 320. Using numer- 
ical models we derive a ^He/H abundance ratio by number of 1.9 x 10^^ with an uncertainty 
of roughly a factor of two. This value is two orders of magnitude higher than the typical ^He 
abundance found in Galactic H ii regions. This detection makes J 320 the second planetary 
nebula known (together with NGC 3242) to have an anomalously high ^He abundance. 
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Fig. 1. — '^He/H abundance as a function of Galactocentric radius where pHe/H] = 
log(^He/H) - log(1.5 X 10~^). The sohd points and crosses are Galactic H ii region ^He/H 
abundance ratios taken from Bania et aL (2002). The crosses correspond to H ii regions 
where no ionization correction could be made and thus we may have underestimated the 
^He/H abundances in these objects. The typical error is shown in the lower right hand 
corner. The triangle plots pHe/H] for the planetary nebula NGC3242 (Balscr et al. 1997). 
The thick solid line indicates the primordial ^He/H abundance predicted by Big Bang nu- 
cleosynthesis (Buries et al. 2001) assuming the WMAP baryon-to-photon ratio (Spergel et 
al. 2003). 
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Fig. 2.— VLA continuum image of NGC6572 at 8.7 GHz. The resolution, 8'.'20 x 7'.'49 at 
a PA of —14°, is indicated in the lower left-hand corner. The intensity scale is in units 
of mJybeam"^ The contour levels are at —3, 3, 25, 100, 500, and 1000 times the 3a 
(0.834 mJy) noise level. 
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Fig. 3.— VLA continuum image of J 320 at 8.7 GHz. The resolution, 7('82 x 7'!51 at a PA of 
— 2 ° , is indicated in the lower left-hand corner. The intensity scale is in units of mJy beam"^. 
The contour levels are at —1, 1, 3, 5, 10, 25, and 50 times the 3(7 (0.264 mJy) noise level. 
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Fig. 4. — H91a; spectrum of NGC6572. The solid histogram is the VLA data toward the 
peak continuum position. The dashed curve shows the NEBULA model assuming Vlsr — 
8kmsec~^. The modeled line at negative velocities is the He 91a transition whereas the 
modeled line at positive velocities is the H 154e transition. The rms over the line-free spectral 
regions is 500/iJy beam~^. 
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Fig. 5. — ^He^ spectrum of NGC6572. The solid histogram is the VLA data toward the 
peak continuum position. The dashed curve shows the NEBULA model assuming Vlsr = 
8kmsec~^. The modeled line at negative velocity is the HlTlr/ transition. The rms over 
the line-free spectral regions is lOOyuJy beam""^. The model '^He'''/H"'' abundance ratio is 
1.0 X 10~^ by number. 
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Fig. 6. — H91a spectra of J 320. The solid histograms are the VLA data and the dashed 
curves are the NEBULA model results assuming Vlsr = — 54.9kmsec~^. Left panel: spec- 
trum toward the peak continuum position for a synthesized HPBW of 8". The rms over 
the line-free spectral regions is 50/iJy beam~^. Right panel: spectrum toward the peak 
continuum position for a synthesized HPBW of 20". The VLA image was convolved with 
a circular 20" beam using the task CONVL. The rms over the line-free spectral regions is 
70 /iJy beam~^. 
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Fig. 7. — '^He'^ spectra of J 320. The solid histograms are the VLA data and the dashed 
curves are the NEBULA model results assuming Vlsr = — 47.9kmsec~^. The VLA data 
were convolved with a circular beam using the task CONVL. The HPBW is shown in the 
top left corner of each panel. The rms values over the line-free spectral regions are 70, 110, 
150, and 190 yuJybeam"^ for HPBW's of 8, 20, 30, and 40", respectively. In the latter three 
cases the estimated S/N is 4. The model ^He'''/H"'' abundance ratio is 1.9 x 10^^ by number. 
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Fig. 8. — Cumulative ^He^ flux density versus radius for J 320. The task IRING calculates 
the integrated flux density in concentric annuli from a specified center location. The solid 
curve is the cumulative flux density for the ^He"*" image using the observed peak continuum 
emission as the center location. The cumulative flux at larger radii drops because of the 
missing short spacings. The dashed lines are the result of a similar analysis using the 
modeled data assuming a halo size of 35" and an electron density of 50 cm^'^. The model 
^He^/H"*" abundance ratio is 1.9 x 10"'^ by number. The missing short spacings do not affect 
the model. The dotted curve shows the result of a similar analysis made for ^He"^ emission- 
free data channels. The rms of the cumulative flux in this noise image within a radius of 20" 
is 43/iJy whereas the peak cumulative ^He"^ flux is 400 /iJy yielding a. S/N for the line of 
~ 9. 
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Fig. 9. — Cumulative ^He"*" flux density versus radius for J 320. The solid curve is that of 
Figures. The dashed lines are the results of a similar analysis using modeled data. The 
models have a halo electron density of 100 cm~^ and different halo sizes from left to right: 
0, 20, 30, 35, and 40". The model results were scaled to match the peak ^He"*" cumulative 
flux density. 
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Fig. 10. — Cumulative continuum flux density versus radius for J 320. The solid curve is 
the cumulative flux density for the continuum image using the observed peak continuum 
emission as the center location. The dashed lines show the result of a similar analysis using 
modeled data. The models have a halo size of 35"and different halo electron densities from 
bottom to top: 0, 10, 25, 50, and 100 cm'^. 
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Table 1. VLA Observational Parameters 



rarameters 


NGC6572 


J 320 


Date 


1 QQ5 Februarv 


1996 Aueust 


Total Time (hr)^ 


24 


23 


Configuration 


D 


D 


R.A. of field center (B1950) 


18''09™40 57^ 


05'^02™48 6"" 


Dec. of field center (B1950) 


06° 50' 25^'5 


10° 38' 25('0 


FWHM of primary beam 


5^4 


5^4 


PWITM of svnthesised beam 


8"20 X 7"49 


7"82 X 7"51 


PA of beam 


-14° 


-2° 


Total bandwidth (MHz) 


6.25 


6.25 


Number of channels 


31 


31 


LSR central velocity (kmsec~^) 


8.0 


-37.9 


Spectral resolution (kHz, kmsec"^) 


195.3 (6.8)'^ 


326.2 (11.4)^= 


Flux density calibrator 


B1328+307/B0134+329 


B1328+307/B0134+329 


Phase cahbrator 


B1730-130/B1749+096 


B0446+112/B0316+413 


Line channel RMS ( //Jy beam"^)*^ 


100 


70 


Continuum RMS ( //Jy beam"^)*^'^ 


280 


90 



^Approximate time not excluding edited data. 
''On-line Hanning smoothed. 
■^Off-line Hanning smoothed. 
'^Using the ^He"^ band. 
^Dominated by dynamic range. 



Table 2. VLA 8.7 GHz Continuum Properties^ 





a 


5 


Peak Flux Density 


Integrated Flux Density 


Angular Size FWHM 


Source 


(1950) 


(1950) 


(Jy beam"''^) 


(Jy) 


(arcsec) 


NGC6572 


18 09 40.54 


06 50 25.4 


0.97 


1.25 


8.4 X 9.4 


J 320 


05 02 48.35 


10 38 21.0 


0.018 


0.024 


8.3 X 9.4 



Based on two-dimensional Gaussian models. 



Table 3. Planetary Nebula Model Parameters 





dinner 


^outer 




T 










Source 


(arcsec) 


(arcsec) 


( km sec" ^) 


(K) 


(cm-3)io3 


(4He+/H+) 


(4He++/H+) 


(3He+/H+)103 


NGC6572 


2.0 


7.9 


15.0 


10300 


22.0 


0.11 


0.00 


1.0 


J 320 


2.0 


7.4 


46.0 


7500 


1.5 


0.10 


0.00 


1.9 


J 320 (halo) 


7.4 


35.0 


25.0 


7500 


0.05 


0.10 


0.00 


1.9 



